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For the steady-state direct cascade of two-dimensional (2d) Navier-Stokes turbulence, we derive 
analytically the probability of strong vorticity fluctuations. When zu is the vorticity coarse-grained 
over a scale R, the probability density function (PDF), Vivo), has an asymptotic behavior lnP ~ 
—zu/zurms at zu 3> zu r ms — [H \n(L / R)] , where H is the enstrophy flux and L is the pumping 
length. Therefore, the pdf has exponential tails and it is self-similar, i.e. it can be presented as a 
function of a single argument, zu/zu rm3 , in distinction from other known direct cascades. 
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I. INTRODUCTION 



Turbulence is a paradigmatic far-from-equilibrium state of matter and the central question of 
physics of turbulence is that of universality: how much one needs to know about an external forcing 
(or initial data for decaying turbulence) to predict the basic features of flow statistics. A related 
question is that of symmetries of the statistics, particularly whether scale invariance appears for 
the scales distant from L, where turbulence is excited [l|, |2j. One distinguishes direct and inverse 
cascades occurring at the scales much smaller or much larger than L, respectively. Data suggest 
that the statistics of inverse cascades are scale invariant [2|-|4|. For example, the probability density 
function (PDF) V of w that is the vorticity u coarse-grained over a scale R is empirically found to be 
a function of a single variable rather than two in 2d inverse energy cascade: V(w) = w~ l f(wR- 2/?J ) 
[7-lo||. Such self-similarity was never observed in direct cascades for whatever small R/L, the 



probability distributions of w change their forms with varying the ratio R/L 0-0. 

One way to explain this profound difference is to argue that fluid motions are slower when the 
scales are larger. As an inverse cascade proceeds, it has an ample time to be effectively averaged over 
the small-scale fluctuations including those of the pumping, whose only memory left is the value of 
flux it generates. On the contrary, small-scale fast fluctuations in a direct cascade stay sensitive to 
the statistics of fluctuations at larger scales nonlinearity then enhances the effect of fluctuations 



>■ . down the cascade so that the small-scale statistics is dominated by rare strong fluctuations. 

One can also explain the difference between the direct and inverse cascades using the Lagrangian 
language. Correlation functions are accumulated along Lagrangian trajectories (for the forced turbu- 

- - ■ lence) or originate from initial data transported along the Lagrangian trajectories (for the decaying 
turbulence). Correlation functions are then proportional to the time the trajectories spend within 
the volume of size L. When the trajectories approach each other back in time, correlations appear at 
larger and larger scales, which corresponds to an inverse cascade. In this case, two-particle behavior 
effectively determines evolution of multi-particle configurations and the second moment exponent 
determines the scaling of higher moments. Just the opposite, direct cascades correspond to trajecto- 
ries separating back in time, one then relates the breakdown of scale-invariance at vanishing viscosity 
to non-uniqueness of explosively separating trajectories in a non-smooth velocity field; exponents 
of higher moments are then related to the law of decay of the fluctuations of the shapes of multi- 
particle configurations which depend on the number of particles jl, Jj)]. Say, for the passive scalar, 



the exponents are independent of the scalar pumping but are dependent on the mixing velocity 
statistics [l3j. Still, one needs to know an infinite number of forcing-related parameters to predict 
the scalar statistics at however small scales. For the direct cascades in Navier-Stokes turbulence and 
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similar nonlinear problems, it is not even known whether the exponents of the velocity moments are 
universal or not. 

Prior knowledge was based on experimental and numerical data, the only analytical results were 
obtained for passive fields in synthetic flows jH, 14-17|. Here, for the first time, the vorticity PDF 



tail is analytically derived from the equation of motion. 

We consider the direct (enstrophy) cascade of 2d turbulence 18l-l20|. which in Lagrangian terms 
is peculiar since it corresponds to an exponential separation of trajectories. Indeed, the physical 
mechanism of the cascade is that pumping-produced vorticity blobs are deformed by the flow into 
thin streaks: stretched in one direction and contracted in another one until viscosity dissipates them, 
like for a passive scalar. An important distinction from the passive scalar is that the vorticity u 
is related to the flow velocity v: u = curl v. Constancy of the enstrophy flux over scales, ((v\ ■ 
Vi + v 2 ■ V 2V1W2) = const, suggests the scaling \v(r) — v(0)\ oc r i.e. spatially smooth velocity. 
In a steady state, the enstrophy dissipation i/|Vu;| 2 must stay finite in the inviscid limit v — > 0. 
The velocity then cannot be perfectly smooth, but the possible singularities are no stronger than 
logarithmic [HI, S3 • If one assumes self-similarity in a sense that the PDF of the coarse-grained 



vorticity is V{w) = w~ l f[m a / ln(L/i?)], then the flux constancy requires a = 3 [18l. l2ll. 123 

There are some consequences of the self-similarity. Say, the enstrophy transfer time through a given 
scale R, determined by the stretching/contraction rate, can be estimated as a turn-over time or an 
inverse vorticity at this scale. On the one hand, that time decreases with the scale as ln~ 1 ^ 3 (L/i?), 
which would suggest that the small-scale statistics is sensitive to the statistics at larger scales. On 
the other hand, the total time of enstrophy transfer from L down to the viscous scale 77 diverges 
oc ln 2 ^ 3 (L/rj) as 77 -> 0. Particle trajectories are then expected to separate exponentially rather than 
explosively and stay unique even in the inviscid limit, that makes self-similarity plausible, according 
to the above Lagrangian arguments. Note that von Neumann (24| and Kraichnan [l8| argued that 
an infinite number of vorticity conservation laws can make the vorticity cascade non-universal, yet 
Falkovich and Lebedev later argued that the fluxes of higher vorticity invariants must be irrelevant 



due to the phenomenon of "distributed pumping" 21| . Recently, self-similarity breakdown was found 



empirically for the vorticity isolines, which are conformal invariant in the inverse cascade, while in 
the direct cascade they are not scale-invariant but multi- fractal with the fractal dimension 3/2 and 
higher dimensions saturating at 1 0,9 ( that 

may be related to strain persistence that leads to 
vorticity organized in long thin streaks). That makes it natural to expect that the bulk vorticity 
statistics is not self-similar as well. 

The present work is devoted to analytical description of a single-time vorticity statistics in the 
steady-state 2d turbulence in the direct cascade. We analytically derive the non-Gaussian tail of the 
PDF of w, that is the vorticity coarse-grained over the scale R, in the direct (enstrophy) cascade. 
We show that the tail is exponential, 

hp(D »~-pHl|p' (L1) 

for a driving force with a finite correlation time. In particular, Eq. (11. ip shows that the vorticity 
PDF is self-similar, i.e. it can be presented as V{vj) = w" 1 f[w z / \n{L/ R)]. Moreover, up to the 
order-unity factor, the tail (II. ip is determined by a single parameter, H, that is the flux of the 
squared vorticity, which also determines the w rms i.e. the bulk of the pdf. To obtain the single-point 
vorticity PDF, the ratio L/R should be substituted by y/Re in Eq. (11. ID . 
The structure of our paper is as follows 
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II. BASIC EQUATIONS 

The incompressible 2d Euler equation can be written for the vorticity uj = d x v y — d y v x : 

^L + {vV)u> = <f>. (2.1) 

Here <fi is curl of the external force / exciting the turbulence: <fi = d x f y — d y f x . The viscous term 
is omitted in f!2.ip . which means that we consider flow variations on scales much larger than the 
viscous scale rj. We shall describe the flow in the Lagrangian reference frame attached to a fluid 
particle placed at the origin, such that v(0) = 0. Then the velocity is expressed via the vorticity as 



-a/3 



4 + 4M-<). (2.2) 



27r If- r' 2 \r 



The pumping is assumed to be a random Gaussian field spatially correlated on the scale L and 
short correlated in time. Then its variance is (0(0, O)0(t, r)) = 5(t)x{r), where x{ r ) rapidly tends to 
zero as r exceeds L. As we shall see below, the processes that contribute to the vorticity PDF tails 
take a long time which allows effective averaging over forcing so that our results are asymptotically 
valid for any forcing with a finite correlation time. 
The statistics of the flow can be examined within the framework of Martin-Siggia-Rose formalism 



25H28| so that all the averages (correlation functions) characterizing the flow are calculated as 



functional integrals, J T>pT>u exp(iX) . . . , with the effective action 



dtd 2 rp{r) d t uj + v'Vuj + - j d 2 r' x(\r — r'\)p(r') 



(2.3) 



Here p is an auxiliary field introduced to put the equation of motion (12.11) into the exponent. Since 
the action (12.31) contains a cubic term originating from the nonlinear term in Eq. (12.11) . one is unable 
to calculate the functional integrals explicitly. Nor one is able to treat the third-order term by a 
perturbation theory, since there is no small parameter in the expansion. In other words, we deal 
with the theory where the coupling is strong. What allows for an analytic description is that we 
consider rare strong fluctuations i.e. describe tails of the vorticity PDF. 

We consider the PDF of w, that is the vorticity u coarse-grained over a scale R from the interval 
of the direct cascade, that is we assume that R is much smaller than the pumping scale L but larger 
than the viscous length r\. A general strategy to find tails of the PDF, V(w), is to calculate the 
corresponding functional integral in the saddle-point approximation utilizing the ratio vj/w rms as 
a large parameter. The way to do that is the so-called instanton formalism adapted for turbulence 
problems 0, io|-[33|. In this way, one looks for an extremum of the action (I2.3p . defined by the 



instanton (extremum) equations SX/Soj = = 5X/5p with appropriate boundary conditions. Both 
the action and the measured quantity u are invariant with respect to rotations and so are instanton 
equations and their boundary conditions. However, axial symmetry turns nonlinear terms in the 
instanton equations into zero killing dynamics. In other words, a "naive instanton" is meaningless. 
The physical reason is quite transparent: there is neither stretching nor contraction for axially 
symmetric flows so that the force can pump the vorticity forever. That means that flow realizations 
that determine a given large value of w must have their axial symmetry broken. We establish below 
that the angle- dependent part of the vorticity realizations remains much smaller than the isotropic 
part during most of the evolution (by virtue of the large parameter w /w rms ). That will allow us to 
integrate over angular degrees of freedom (in the Gaussian approximation) and obtain a renormalized 
action for the zero harmonic Uq. Moreover, we show that only the second angular harmonic provides 
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for the relevant renormalization by virtue of the large parameter ln(L/r). We then find the new 
(effectively axially symmetric) instanton that corresponds to the renormalized action and gives the 
tail of the coarse-grained vorticity PDF, Vivo). 



III. SEPARATION OF HARMONICS 



We use polar coordinates, x = r cos<£>, y = rsin<£>, and expand the fields u and p over the angular 
harmonics: 

u(t, r) = }^cu m (t, r) expjirrup), 2irp(t,r) = ^Pmjt, r) expjirmp). (3.1) 

Then the effective action ( 12. 3 p splits into a number of terms Z = l +J2 m>0 (T m +T_ m +Z mm +I im )+I 3 , 
where Z contains only the zero harmonics ooo^Po- The last term Z 3 is a sum of the third order terms 
containing harmonics with m ^ 0, it is neglected in what follows, which is justified below. The terms 
quadratic in p m ,u m are written as follows 

Z = J dt dr r p d t uJo + - J dt dr r dr' r' xo( r ^ r ')Po( r )po( r ') , (3.2) 

l m = J dt drrp_ m {d t u m + v imu rn /r + d r u v rm } , (3.3) 

l im = - J dt dr r d r p {v r , m U- m + v r ,- m u m ) , (3.4) 

Z mm = i J dtdrr dr' r' Xm(r,r')p-m(r)p m (r') , (3.5) 

Xm(r,r') = J^ eimV X + r' 2 -2rr>cosip^) . (3.6) 

Here v in Eq. ( I3.3P is related to ujq via the equation ojq = v /r + dv /dr. Our goal now is to derive 
an effective action for the zero harmonic, Z + AZ by integrating over all the other harmonics, 

exp (iZ + iAl) = [ T>uj ±m Vp± m exp (zZ) . (3.7) 

m>0 

The integration is Gaussian if to neglect third order terms in p m and cu m , as explained above. Then 
AZ = Yl AZ m where 



e iA± m = j T>u± m Vp± m e llm+tl - m+llmm+ll > m . (3.8) 

If I im = then the expression ( 13.81) is the normalization integral that is equal to unity due to 
causality j28|. Therefore one can write 

OO - n — l 

AJ = EEV ((Xm)n)c ' (3 - 9) 

where the angular brackets mean integration over u± m and p± m with the weight exp(zZ m + zZ_ m + 
iX mm ) and the subscript c means an irreducible average (represented by connected diagrams). 

Consistently considering small fluctuations (as in neglecting Z 3 ) we take only the term with n — l 
in Eq. ( I3.9p . We shall justify it later by observing that the angular part remains small during the 



build-up of the strong fluctuation that we consider. Therefore, the main object, that contributes to 
the n — 1 term and need to be examined, is the pair correlation function 



(uJ m (t, r)uj_ m (t', r')) = F m (t, t'; r, r'). 
The simultaneous pair correlation function satisfies the equation 

vo(ri) v (r 2 ) 



d t F m (t,t,r 1 ,r 2 ) + im 



ri 



F m (r 1 ,r 2 ) 



-d r u (ri)- I drr 2% 



\m+l\ 



U 



\m-l\ 



r — n 



-F m (r, r 2 ) 



+d r uj (r 2 )- I drr 



2 U 2 



\m+l\ 



\m— 1 



\T — Tn 



-F m {r u r) = Xm(ri,r 2 ) 



(3.10) 



(3.11) 



where U\ y2 = min{r/rx 2 , T 12 /r}. One should treat separately the first angular harmonic, with 
m = ±1. 



A. Logarithmic approximation 



Let us pass to the logarithmic variable £ = ln(r/L), where L is the pumping length. We are 
interested in small scales, r < L where |£| ^> 1. We consider only the leading contributions in terms 
of large |£|. In this case, only the terms with m = 2 are relevant in Eq. (I3.9P since the integration 
in the expressions (I2.2f3.lip is logarithmic only for them, as has been noticed already in 21] , the 



feature is likely related to peculiarity of elliptic vortices in straining flows |29j. Other harmonics 
behave as r m ~ 2 for m > 2 i.e. are suppressed exponentially in terms of the large logarithm £. In the 
logarithmic variables, Xm(£i> £2) for m 7^ are nonzero only if both l^l, |£ 2 | ~ 1, since the integral 
(13.6j) is zero for r\ +r 2 < L and decays as Lf ^r\r 2 = exp[— (£1 + £ 2 )/2] for ri,r 2 —> 00. That means 
that one can approximate x m ((i)6) ~ H m 5(^i)5(^ 2 ) for m/fl. The zeroth harmonics can be taken 
as Xo(£i)£2) = -^#( — £i)#( — £2) where 9 is the step function. 

Next, we pass to the field #(£) = r 2 po(r). Then the bare action ( 13. 2 p for the zeroth harmonics is 
rewritten as 



X = J dtdiqd t UJ Q + l - J dtdii d£ 2 Xo(&, 6)<?(6)<?(6) . 



(3.12) 



The correction AX (13. 9 j) in the main approximation (taking into account only the term with n — 1, 
m = 2) can be written as 



AX^i Jdtdtq(t,Z) Jd( [F 2 (tM,t)-F 2 {tM,0] ■ 



(3.13) 



Here the function F 2 has to be treated as a functional of loq to be extracted from the equation (13. lip 
for m = 2. 

In terms of the logarithmic variable £ the equation (13.111) is rewritten as 



$^(£1,6) + 2zK&) - w(£ 2 )]F 2 (£ 1; £ 2 ) -x 2 (6,6 



z'<%j (£ 2 



00 



^2(6,0+ / ^e 4(e - 6) F 2 (e 1 ,e) 



(3.14) 



/ ci£F 2 (£,£ 2 ) + / d£e 4 ^F 2 (£,£ 
J— 00 



where we denoted w(£) = vq/t = ff d£' exp[2(£' — £)]uo(£'). In the main logarithmic approximation 
we get from Eq. (13.141) 

dtFzfa, £ 2 ) + ^ (t, 6)^(6,6) + 0^0(^6) / d£F 2 (U2) 

1 hx 



-iu (t, 6)^2(6, 6) - ^ofrk) / ^^2(6,0 = x 2 (6,6) 



(3.15) 



Putting in Eq. ( 13. 15ft £1 = £2 = £ an d substituting the result into the expression ( 13. 13ft one obtains 

/ dtd£-^-c> t F 2 (£,£), (3-16) 
the term with \2 is neglected since it is small in the interval of the direct cascade. 



B. Eigen functions 

The equation ( 13. 15j) can be rewritten in the form 

^2(6, £2) +i0 1 F 2 -i6 2 F 2 = X2(£i,£ 2 ) , 

where 



o/(£) = J d( 

Let us introduce eigen functions of the operator O 

iaw„(0 



^(£-C) + 2^o^(C-£) 



/(C). 



C^ A = w (£)^a(£) + - 



2 <9£ 



dC^A(C) = A^(£). 



(3.17) 



(3.18) 



(3.19) 



We assume (in accordance with the answer obtained) that Wo(£) monotonically diminishes from some 
value s at £ = —00 to zero at £ = +00. Then a set of the eigen functions of the operator O can be 
written as 



(fx = 0{u Q - X)2d^u - X = d 5 u , 



(3.20) 



where 9 is the step function and < A < s. The functions ( I3.20p are the right eigenfunctions of the 
operator O. Analogously, one can define the left eigenfunctions: 



</v(£) = 7T liTC k Re ^ ~ + ie \ 3/2 

where, again, < /x < s. The functions (13 . 2 1 f) satisfy the equation 



d(d c uj (pf,(C) = /i0 M (£). 



(3.2i; 



(3.22) 



The only thing which is important for what follows is the orthogonality and normality of the right 
and left eigenfunctions. The factors in ( I3.20|l3.2ip are chosen to ensure the normalization condition 



that can be checked directly: 



^A(0^(0 = -7^1imRe ' 



27r e^o J x y/uj - \ (jj, — u + ie) 3 / 2 
= -Vs - AlimRe(A - fi - z'e) _1 (/i + ie- s)~ 1/2 = 5U - A), (3.23) 
provided < /i, A < s. One can check completeness of the set (I3.20|3.2ip : 

d\MOMO = ^-0- (3-24) 



The check is reduced to an integral analogous to one (13.231) . 
Let us expand F 2 over the eigenfunctions (13.201) . 

F 2 (t,(i,(2) = J d\ 1 d\ 2 $(t,\ 1 ,\ 2 )tp Xl (( 1 )tp X2 (( 2 ), (3.25) 

where 

*(Ai,A 2 ) = J ^6^(6)^(6)^(6,6)- (3.26) 
Now we take into account that u>q — > as £ — > +00 and obtain: 

[d t + i(iix- /i 2 )}^(t,ii l ,fi 2 ) + J c/Ai$(Ai,/i 2 )J(/ii,Ai) 

+ J rfA 2 $(/i l! A 2 )J(/i 2 ,A 2 ) = J rf6^2 w (6)^ 2 (6)X2(6,6), (3.27) 

where 

j(/i,A) = [ dCMOdtMC) = - [ d(d t MC)MC)- (3.28) 



The equation ( 13. 27ft is equivalent to Eq. ( 13.171) . 
Substituting into the definition (13.281) the explicit expressions (13. 20113. 21 j) we get 

JM ^^U. x) r X ,3, 9) 

^/i 2 [ Jx V( w o - X)(jjl-u ) J 

Here we introduced the designation 

d t uj (t,Z) =V(f,Wo). (3.30) 
Note that it follows from the definition (13.301) that 

* =^r = dt&7 = dMd ^ ) - (3 - 31) 

Performing the substitution uq = A + (/x — A)x we get from Eq. (I3.29j) 

J(lM, A) = A [*(,, - A)^(A)] + ItfOi - A)V'(A) 

+-0(ijl - A) / ^=V" [A + 0* - A)x] . (3.32) 

We see that J is the sum of singular terms and the term, which can be written as a regular expansion 
over fi — A. 
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C. Initial Condition 



As we will see, the instantonic solution for ooo diminishes back in time. Therefore the saddle- 
point approximation ceases to be correct at a time t = i* where uq is of the order of typical (rms) 
fluctuation and one should pose the initial condition for the instanton at t — t*. We assume that the 
initial fluctuation (at t = £*) uq is some slow (logarithmic) function of the distances in the region 
|£| < ln(L/r) so that d^\nu ~ Next, we assume that at t = t* 

(w(n)w(r 2 )) ~ (F 7 ) 2/3 , (3.33) 



7 = In — — = - In Ur± + r\ — 2r\r 2 cos </?) / L 2 ] 

L 2 

2rir 2 cos ip 
r\ + r 2 



iln 


"r 2 + r|" 


1 . 
+ - In 


L 2 


2 


2 



(3.34) 



It is important to stress that, strictly speaking, we cannot derive the second moment (13.331) from 
the equation of motion. That choice is consistent with the flux relation and, as we show below, is 
self-consistent with the higher moments described by the PDF tail to be derived. 
We now derive the second angular harmonic of the pair correlation function: 

^(4,6,6)= / ■^-exp(-2i<p){u(r 1 )u(r 2 )). (3.35) 
Jo 27r 

If |6|, |6l ^ 1 then a power-like function f(ip) can be expanded in the ratio of the two last terms 
in Eq. (13.340 . Then one obtains 

r?exp(-2^)/( 7 ) = 4-ix r 'lf\^>^ (3.36) 
J 2vr FK TOU/ 4rf 7o r 2 /r 2 if n < r 2 , v ; 

where 70 = (1/2) ln[(r 2 + r 2 )/L 2 } and we used the relation 

f 2n 7T / 

/ dx cos(2x) ln(l — a cosx) = — ( a 2 — 2 + 2yl — a 2 ) . 
Jo a v / 

We conclude that the function F 2 (t*) exponentially in 6 — £ 2 tends to zero as |6 — 61 = I m (ri/L) — 
ln(r 2 /L)| tends to infinity. Therefore it can be estimated as F 2 (Q ~ i7 2/3 |6|- 1/3 5(6 - 6)- 
Now we can analyze the initial value of the function $ 



$(t„Ai,A 2 ) = J ^6 ^(6)^(6)^2(4, 6, 6)- (3-37) 
Substituting here the above estimate, one obtains 

$(t„ Ax, A 2 ) ~ H 2 ' 3 f dC\C\- 1/3 K(0K(0- (3-38) 
«/ —00 

Again, the approximation (I3.38j) implies Ai, A 2 3> H 1 / 3 . What is most important is that the integrand 
of ( I3.38P has a third-order pole when Ai — > A 2 — > ujq. Integration over £ results in the second-order 
pole: $ (Ai, A 2 ) oc (Ai — A 2 )~ 2 at Ai — > A 2 . Therefore 

(HE) 2 / 3 

*{U, Ai, A 2 ) ~ i-^— ^(Ax - A 2 ), (3.39) 
v AiA 2 



where £ is determined by the condition wo(£*,£) — Ai ~ A 2 . 
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D. Adiabatic approximation 



The instantonic field u;o(i, £) describes an optimal fluctuation that starts from the rms level at 
the time t* and grows to a prescribed large value of w. After we find below the instanton solution 
a>o(t, £)> we see that its form changes slow on its own rotation timescale uJq 1 - In particular, that 
means that the eigenf unctions ( I3.20p change slow too, so that one can neglect non-singular term in 
(Frm ). Then 



the expression 



JGu, A) -+ A [8(im - A)^(A)] + ^ - A)^'(A) 



(3.40) 



Substituting the expression (13.401) into Eq. (13.271) and omitting the right-hand side we get 

[d t + i(jJLl - /i 2 )] ®(t, /J.1, /U 2 ) + "0(^1 



— $ + 4>Mq^$ + g + = . (3.41) 



The equation (13.41 j) can be solved by the method of characteristics. The equation for the character- 
istics reads 



~dt 



ip{t,fi) = d t u (t,£) 



(3.42) 



where £ is a function of \i to be extracted from the condition \i — ouq- An obvious solution of the 
equation (13.421) is 

^(t)=w (t,0, (3-43) 

where £ plays the role of the marker of the characteristic. A solution of the equation (13.411) is written 

as 



x exp | 
Using Eq. (13 . 3 1 1) we get 



$(t,/ii,/i 2 ) = $[t*,/ii(t*),/i 2 (t*)] 

3 3 
-z>i(s) + i/i 2 (s) - -^'[s,/ii(s)] - ~^'[s,// 2 (s)] 



(3.44) 



a c Wo(t*,Ci) 9 c w (t*,C2) 



_<9 ? wo(t,Ci) <9 ? cu (t,C 2 



3/2 



x exp 



ds [-i/i x (s) + «/U 2 (s)] 



(3.45) 



where the variables d, C2 have to be extracted from the relations fii = co (t, Ci), /i 2 = w (i, C 2 )- 

Now we can estimate a role of the regular contribution omitted in Eq. (I3.40p . see Eq. ( I3.32p . 
Substituting the expression (I3.45P into Eq. (I3.27P we then conclude that an integration over Ai or 
over A 2 in the omitted terms is determined mainly by the oscillating factor in Eq. (13.451) . Then 
J d\\ — y (t — t*)" 1 (and the same for the integration over A 2 ) and therefore the omitted terms in the 
equation (13 . 2 7f) can be estimated as 

1 

y (/i!)$0ui,/i 2 ) . 



t - 1 



Comparing the term with Eq. (13.41 p . we conclude, that the omitted terms have an additional small 
factor [{t — t^uo]' 1 and therefore the approximation leading to Eq. (13.411) is correct. 
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Let us calculate F 2 (£i, £2) using the expression (I3.45P : 

F,( t , 6 .6)=/W^-^- , ' ll y;-^-^ (3-46) 
xfi n (f P \F> n (f P \ P^Q^Ci) d^ofo., C2) 1 ^ 

x$[t»,o;o(t*,Ci) > Wo(t*,C2)]exp [-iw (s, Ci) + «*>o(s, C2)] 

Recall, that £1 and £2 have to be extracted from the relations /ii = cj (t, Ci); ^2 = ^0(^2)- After 
substituting (I3.39P into (13. 46 p . in the integral over /i 2 we keep only the pole term at ji 2 —> 
£2 — > Cij differentiating the exponent gives the t — factor because of the slowness of the instanton. 
Passing then from the integration over fii to ones over ( we obtain 



F 2 (t, e, ~ H 2 '\t - U)[d(U> (t, 0? f d( ( 2/3 d c In u> {U, () 



9[uo{t,Z)-a)o{t,0] 



wo(*,0 -w (t,C) 
~ (t - t„) (^ )^ 2/3 r 1/3 ln(w /^o) . (3.47) 

Here we used 9^1no;o(t*,C) — 1/Cj ^ is equivalent to assuming that the initial u; (£) is a power- 
like function of the logarithm £. The logarithmic divergence in (I3.47P is cut off due to a finite 
(order-unity) width of F 2 (t*, £1, £ 2 ) over £1 — £ 2 - 
Now we turn to the pumping contribution: 

pt-U 

$(n,\,t)= dr exp[-i(p- A)(t-t* -r)]S(//,A), (3.48) 
•/ 

S(//,A) = y ^i^2 M (ei)0A(6)X2(ei,6). (3.49) 

The function X2(£i ; £2) is nonzero provided £i,£ 2 ~ 1- Since the integration over ^ and £2 in Eq. 
(I3.49P smears the singularities in then S a smooth function of \x and A. Next, a characteristic ujq 
in the integral lEElgj) is if 1/3 . Therefore S oc /i _3/2 at /i > iJ 1/3 and H oc A _3/2 at A > i/ 1 / 3 . Thus 
the characteristic [i and A in the integral (13.251) are less or of the order of H 1 ^ 3 . Then we obtain 



Fi ™ V{u) w (W /■ d dX r* dr 

y Jo 



exp H(/i - A)(t - t, - r)] S(/i, A) ~ # 1/3 1^1!. (3.50) 

Since we shall obtain a slow instanton with if 1 / 3 ^ 3> 1, the contribution (I3.47P is larger than (I3.50p . 
That means that the pumping-produced anisotropic fluctuations give lesser contribution than defor- 
mation of an initial fluctuation. The consequence is that the tail of the vorticity PDF is insensitive 
to the form of the pumping correlation function and is determined solely by its zeroth moment i.e. 
the vorticity flux. That means universality of the statistics of strong vorticity fluctuations. 

The estimation (13.471) determines the main contribution to F 2 (^i, £2) where £1 ~ £ 2 . We checked 
this time dependence of F 2 by solving numerically Eq. (13.151) using different time- independent uj . 
At the beginning, we have chosen F 2 being determined by a typical fluctuation. Then we checked 
that F 2 linearly grows as time runs. It is interesting to note that the antisymmetric in £x, £ 2 part in 
F 2 saturates. The behavior is in accordance with Eq. (I3.15p . 
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IV. INSTANTON 



Substituting the expression (I3.47P into Eq. ( 13. 1 6[) one obtains 

AX~ -H 2 '* J dtdiqC l, \d^Y 1 dA{t-U)d^\n{uj () /d^)}. (4.1) 
Collecting the expressions (13. 12114. ip we get finally the effective action 

J cff = J dtd^q {«9 iWo -c# 2 / 3 C 1/3 (^^ j dtdtiid&q{{ii)q{&) 

where c ~ 1 and we substituted Xo(£i>£2) — H6{—^i)6{—^2)- By rescaling q, u and H we can put 
c-> 1. 

X c ff = j dtdiq {d t u - H 2 / 3 r 1/3 (d^ )' l d t [(t - Qd^ H^o/d^ )}} 

\H f dtd^d&q^q^), (4.2) 



2 

We are interested in the PDF of w, that is u coarse-grained over the scale R. In our terms, it can 
be written as 

2 r R 

W= B 2 J drruj °~ e>oMR/L)), ( 4 - 3 ) 

because of the logarithmic character of the uq dependence on r. Thus, we should fix uq[\tli(R/L)} = w 
at the observation time. Since we consider steady-state turbulence, the moment of measurement is 
arbitrary, we choose it to be t = 0. Then PDF V{w) can be calculated as the path integral 

V — Jvloo Vq exp (il cS ) , (4.4) 

taken at the condition u; [0, ln(r/L)] = w and for the fields u , q defined at negative times t M. The 
last property is explained by causality: the values of the fields at positive times cannot influence 
the PDF V{w). Since w is fixed at a single point at t = the field q at t = satisfies q(Q oc 
8[£ — \n(R/L)] that reflects the measuring procedure. 

In the saddle-point approximation, we put \nV ~ i1^ r where X|ff tr is the extremum value of the 
effective action. The extremum conditions 5X eS /5u = = 5Z eS /5q give the so-called instanton 
equations: 

dtooo = H 2 ' z C l/ \d c uj )- l d t [(t - Qd^uo \n(ooo/d c uj )} + HQ(t) , (4.5) 
d t q + H 2 / 3 d c {HtUo/d^oX-^id^oy'dtKt - Qq}} = , (4.6) 

where Q(t) = —i f d(q((,t). In deriving Eq. (14.61) we exploited large value of the logarithm 
ln(uo/d^uJo) ~ In |C| s° that in the main order we only account for the terms in the equations 
that contain the logarithm. Apart from the logarithm, the correction (14. ip depends only on the vor- 
ticity spatial derivative d^UQ. As a result, the variation with respect to vorticity gives the equation 
(14. 6p . which has the form of a continuity equation, so that dQ/dt = in the main order. We see from 
Eq. (14. 5 p that the first term in the right-hand side is negative at £ < that is the correction (14.11) 
describes decrease of the vorticity due to deformation of the circular vortex by elliptic perturbations. 
Substituting the relation (14. 5p into the expression (14.21) . one finds 

\nV w ~ [dtQ 2 ^ -^-Q 2 \U\. (4.7) 



2 / 2 
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Since Q is t-independent, one readily obtains from Eq. (14 . 5 [) that cuo grows linearly with time, 
^o{t, = ' (t — t*). Then we find for the factor 8 



H 2 ' 3 8 

Replacing here the ratio B/d^B by £ one obtains 

#2/3 



U 



£1/3 



ln\Z\ + HQ 



{t-Q, 



w 



#2/3 



[ln{L/R)}V 3 



In ln(L/r) + HQ 



Q 2 w 



Then one obtains from Eq. (14. 7p 

~ ~2{Q -2H- 1 /z[\ n (L/R))- l /nnln(L/R)} 
Optimizing the expression over Q one gets 

\nP ~ -4iJ- 1/3 [ln(L/i?)]~ 1/3 lnln(L/r)u; 

#2/3 



TO 



[ln^/ii)] 1 ^ 



lnln(L/r)|t„ 



(4.8) 

(4.9) 
(4.10) 

(4.11) 

(4.12) 
(4.13) 



The value of u)o(t*) does not matter with logarithmic accuracy as long it is much smaller than w. 
The expression (I4.12p leads to the final answer (11.11) where we omitted the slow factor ln[ln(L/i?)]. 

We can use the instanton solution found to check the validity of all the assumptions made in the 
derivation of the effective action. Remind that we consider the case ln(L/i?) ^> 1. The applicability 
condition of the saddle-point approximation is 



\w 3 \ ^>H\n(L/R) 



(4.14) 



The fluctuations on the background of our instanton are indeed small: using the instanton solution 
we estimate F 2 — t^co ^~^ 3 ujQ/^ luq as was assumed. That justifies neglecting £ 3 and n > 1 
terms in (13. 9p . Let us now estimate F. m for m > 2 and compare it with F 2 . Zero mode of (13.111) 
must allow cancelation of integral and non-integral terms which is possible only when the integrals 
in (13. lip are logarithmic. That requires F m oc r m ~ 2 oc exp[(2 — m)£], therefore those terms are 
exponentially suppressed comparing to F 2 . The instanton duration time is such that \wt*\ ^> 1 so 
that our instanton is "slow", that is indeed u}o(t) changes slowly comparing to itself. 



V. DISCUSSION 

It is illuminating to compare vorticity statistics in the direct 2c? cascade with the statistics of the 
passive scalar in a spatially smooth random flow 12|, ll8N20|. For a passive scalar 9 coarse-grained 
over a scale R less than the pumping length L one can get the asymptotic behavior of the single-point 
PDF in a smooth random flow by the following simple reasoning. Large values of 9 are achieved when 
there is no stretching for a time which is much longer than the mean stretching time A -1 ln(L/R), 
where A is the Lyapunov exponent. During that time, the passive scalar is pumped by a random 
forcing, i.e. it has Gaussian statistics with the linearly growing variance: 



V{0) 



dt Q(t) exp(-# 2 / Pt) 



(5.1) 
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where Q(t) is the probability of no stretching during time t. Stretching is correlated on the velocity 
timescale Aq , which is independent of 9. For every stretching event, the scalar blob is stretched by 
order e and we ask for the probability that there were less than the number In Pe such events during 
t. For t ^> Aq , this is the probability of the Poisson process In Q(t) ~ — \ t + 0\\ia (L/R)}. Doing a 
saddle-point integration over t we obtain the exponential tail (first suggested in |34j, |35[ and derived 



by the instanton formalism in [3JJ, l36|): 



HV{9)\ ~ -9^/\jP + 0[\n{L/R)\. 



(5.2) 



In the work 21 



we established some features of the direct vorticity cascade using an analogy 
between the vorticity u and passive scalar 9. Developing that analogy further, one can propose an 
interpretation of the tail (11.11) . For the vorticity cascade, we may use similar reasoning as for the 
passive scalar with the knowledge added from [2l| that the stretching correlation time is the mean 
total stretching time ln 2 ^ 3 (L/i?)/if 1//3 from the scale R to L. That gives 



~ J dt Q(t,w)exp(-w 2 /Ht) ~ J dt exp -w 2 /Ht -tH 1/3 \n~ 2/3 (L/R) 



(5.3) 



The saddle-point integration shows that the main contribution comes from t ~ w ln 1/3 (L/R)H^ 3 
in agreement with (14.131) . and the result lnf'P^)] ~ — \u\/\H ln(L/i?)| 1//3 reproduces the dependence 
of (II. ip . We see that vorticity is indeed like passive scalar: the stronger the fluctuation the longer 
it lives which leads to a sub-Gaussian PDF tail. Such tails were observed in numerical simulations 
For quantities (like velocity) whose statistics is determined by fast events, their PDF have 



40 



tails steeper than Gaussian [4JJ. 

In a finite box, coherent vortices may appear due to an inverse cascade [42 -44j. The vortices have a 
well-defined spatial profile of the average velocity field that is explained by an interplay of the average 
profile time derivative (or friction) and an effective pumping related to long-correlated fluctuations 
[451 ] . An interesting question that is a subject of future investigations concerns an influence of the 
coherent vortices on the enstrophy cascade. One can also think about extension of our instantonic 
approach to the inverse (energy) cascade of the 2d turbulence. 
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